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\^ . Abstract 

o ■ 

Q\ I In this article, we study the 5-state Potts random matrix models extended to 

0^ ■ branched polymers, by the equations of motion method. We obtain a set of loop 

equations valid for any arbitrary value of q. We show that, for q = 2 — 2 cos -vr 

{I, r mutually prime integers with I < r ), the resolvent satisfies an algebraic 

equation of degree 2r — 1 if / + r is odd and r — 1 if / + r is even. This generalizes 

,^ , the presently-known cases of q = 1, 2, 3. We then derive for any < g < 4 the 

Potts-g critical exponents and string susceptibility. 
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1 Introduction 

Random matrix models |1| have proven to be a powerful mathematical tool for the 
study of statistical physics systems on a fluctuating two dimensional lattice [0]. In 
particular, the Potts model [Q on a random lattice [Q, which was first partially solved 
by Daul 0, has received a recent renewed interest [§, |^ due to new approaches to the 
problem. It is a g-matrices model where all the matrices are coupled to each other, 
which prevents one from using the formula p, ^ to integrate out the relative angles 
between the matrices (they no longer are independent variables) and deal with the 
eigenvalues only. In this paper we use the equations of motion method which does 
not involve integration over angular variables. We obtain non-trivial loop equations 
relating even moments to odd moments of a single matrix Mj, and we show how to 
extend these relations to the case of Potts-g plus branched polymers (gluing of surfaces). 
Such relations (which could not be obtained by previous methods |^, ^) are needed to 
apply the renormalization group method |10[ to Potts-g models with added branching 
interactions, which is hoped to provide an understanding of the c = 1 (g = 4) barrier. 
We obtain an 0{n)-like equation, the solution of which is known [|ll|] and involves 
elliptical functions [Q. What was the resolvent in the 0{n) model, however, is now 
(up to some transformations) the functional inverse of the Potts-g resolvent, with n 
replaced by 2 — g. When g = 2 — 2 cos(z/7r) with u rational, the general elliptic solution 
degenerates into an algebraic function. It has already been observed 0, ^ that for the 
particular cases of g = 1, 2 or 3, the resolvent obeys an algebraic equation of degree 2, 
3, 5 respectively. In this article, we will derive from the value of ly, for any "rational 
g", the degree of the algebraic equation obeyed by the Potts-g resolvent. We will also 
derive the Potts-g critical exponent for general values of g, which agrees with Daul's 
expression. 

2 The Potts-g matrix model 

The Potts-g model with branching interactions (which appear if one wants to apply 
the renormalization group method) is defined by the partition function: 

Z = jdM,. . .rfM,e-^'S»^*^^'+'^(^*'^*^-**^ E.^.A..M,) ^2.1) 

with Mi hermitian matrices N x N. 

The partial derivatives of ip with respect to ti Mf /{2N) and tr J^j^i^iMj/N are 
U and c/2 respectively, and their expectation values (which are numbers) U and c. 
When U = 1 and c is a constant, this reduces to the ordinary Potts-g model. 

We will also define the following functions: 

Wiz)= ^(tr 



NX""- z-Mi 



Wiz)= H^^J^M, 



-^(^' ^ ) ~ AT \ ^^ z-Mi z'~Mj 



(2.2) 

which do not depend on the indices provided i ^ j ^ k. The function F{z, z') and 
F[z^ z') are thus symmetric: 

F{z, z') = F{z\ z) , F{z, z') = F{z\ z) (2.3) 

We will also define: 

f^z) = W{z)-gz^ + {c-U)z (2.4) 

Recall that U and c can be general functions of the numbers {ii Mf) and (tr MjMj). 
The moments tk of the resolvent W{z) are defined by the large z expansion: 

Wiz) ~ - + ^r + . . . + ——- + . . . when z ^ oo 
z z^ z''+^ 

and we define: 

c-U 
u = 



3 Equations of motion 

We are going to work in the large A^ (planar) limit, where we have the factorization 
property: (tiA tiB) = (try4)(tri?) |[T^. The following changes of variables in Eq. ( |2.1| ) 
then give the following equations: 

• 5Mi 



z-Kh ' 
.2t 



g{z'W{z) -z-ti) + U{zW{z) - 1) + c{q - l)W{z) = W\z) (3.1) 

5M2 = \ 



1 1 _| 1 1 



z-Mi z'-M2 z'~M2 z-Mi 

g [z''^F{z, z') - z'W{z) - W{z)) 
+ U{z'F{z,z')-W{z)) 

+ c{zF{z,z')-W{z')) 

+ c{q-2)F{z,z') =W{z')F{z,z') 

(3.2) 



Substracting Eq. ( |3.2| ) with {z <-> z'), and using Eq. ( |2.3| ) we can get rid of F in 
Eq. (P3): 



{f{z) - f{z'))F{z, z') = g (^z'W{z) - zW{z') + W{z) - W{z') - uW{z) + uW{z')) 

(3.3) 
In particular, if we choose z' such that f{z') = f{z) and z' ^ z, then we have: 

{z' - u)W{z) -{z- u)W{z) + W{z) - W{z') = (3.4) 



Eq. ( |3.1D allows us to eliminate W, and we then have an equation involving only 
W or equivalently /: 

/ /(^') = /W 

1 ^{z + z' -u- q^)f{z) = {z + z'f - qzz' - (2 - q)u{z + z') + (1 - q)u'^ - i 

\ c g c 

(3.5) 
This seemingly difficult non-local equation is sufficient to compute f{z) . 

Let us ffist study it perturbatively. When z —^ oo, we have f{z) ~ —gz^, thus 
z' ~ —z. If we expand z'{z) in powers of z by solving f{z) = f{z') perturbatively, 
then insert this expansion into the second equation of Eq. (|3.5| ), we obtain a set of 
equations of motion, the ffist of which are: 

gt2 + {cq - gu)ti = 
g% + g{cq-2gu%- gu{cq-2gu)t2- g{gu^ + 2)ti + c{l-q) +gu = 

(3.6) 

with tk = (tr Mf /N). Let us recall that, in the general case, u and c are functions of 
^2 and ti^i = tr MiMj/N. However, we can compute ti^i thanks to Eq. ( ^.Ij ) in function 
of the tfc's: 

c(g - l)ti,i + gt3 + Uh- 1 = 

Finally, our equations, contrary to ordinary equations of motion, allow us to relate 
even traces of a given matrix Mj to odd traces of the same matrix . 

Only even traces appear as the higher order traces in these equations. Let us explain 
why. If we write the expansion of f{z) as: f{z) = —gz{z — u) + 1/z + Y^'^ ti^i/z\ the 
z~^ coefficient of the f{z') — f\z) = equation reads, as z' = —z + u + . . .: 

ii-iy - 1) ti_i + i-iy-' {t - 1) U.2U + . . . = 

Thus, our loop equations give us a relation between the even and odd parts of f{z), 
and f{z) can be completely determined by the requirement that it has only one cut in 
the physical sheet (i.e. f{—z) is regular along this cut). 

Finally, let us stress that these equations of motion, which are valid for branched 
polymers, are a precious tool whenever one wants to study Potts-g models by the 



renormalization group method ||T0 . 



4 Correspondence with the 0(n) model 



Let us now see how to deal with a non-local equation such as Eq. 

The function z'{z) defined above maps one solution of the equation f{z) = y on 
another. It is involutive in the sense of multivaluated functions. We have: 

z'iz'iz)) = z 

Let zq be a fixed point: 

z'izo) = Zo 



and /o = f{zo). 
Then let us set 

C = V/o - / 

and consider z as a function of (. Then we have z' = z{—() and z{() is regular at 
C = 0. Let us set 

^(C) = ^(0 + -J^iC - /o - (2 - q)cu) 
cA — q 

Eq. ( p.5| ) rewritten in term of uj{C) is now an 0(n)-like quadratic equation: 

.;2(C) + ^^(-0 + (2 - q)uj{Ou{-0 = R{C) (4.1) 

where the right-hand side of the equation is an even polynomial of degree 4. The 
similarity between Potts-g and the 0{n) model had already been noted P, H, but it 
had always been said to be unphysical. We shall show here how one can relate the 
results for the 0{n) model to those for Potts-g. 

Eq. ([4.1| ) can be solved exactly, as in the case of the 0{n) model [0. Here, 



we will assume for simplicity that uj{Q has only one physical cut [a, 6] with ab > 
0. Let us denote q = 2 — 2cos(z/7r), with < z/ < 1. Then we have, by writing 
R{C + iO) -R{C-iO) =0: 

{uj{C + iO) -uj{C-iO)){uj{C + iO)+uj{C-iO) + 2cos{u'K)uj{-C)) = for a < C < & 

Thus we have the linear equation: 

uj{C + iO) + uj{C -i0) + 2 cos{u7i)uj{-C) = (4.2) 

The general solution for uj{() is known and can be expressed with elliptical functions. 
It degenerates, in the rational case (i.e. when u is rational), into the solution of an 
algebraic equation. However, uj{() is not the resolvent of the model as it is in the 
0{n) model. Indeed, it is rather the functional inverse of the resolvent for the Potts-g 
model, up to some transformations. Thus, the phase diagrams and critical exponents 
of Potts-g, as expected, are not the same as for the 0{n) model. 

5 Rational case 

Here we will assume that z/ is rational: u = l/r where / and r are relatively prime 
integers. We first recall how to obtain Eq. (|5.4|) , which is an algebraic equation for 

^(0 O- ^ ^ 

If we denote a;+(C) = e^~ uj{Q + e 2~a;(— ^) and oj-{Q = uj+{—() , then Eq. ([4.1|) 
reads 

cu+i0u^4C) = R{C) (5.1) 

and Eq. ( |4.2| ) becomes: 



uj+{( + iO) = -e'"^ uj_{( - iO) u;_{C + iO) = -e-'"^ LJ+{C - iO) (5.2) 



If 0(C) is defined by: 

uj+ = Vr e'^'t' —^ then we have u. = VR e''^^ —^ 



and 



^(C) 



Rcos{(f)) 



smz/vr 



Eq. (|]2D shows that 



^(c) = ^K + (-ir^^:) (5.3) 

has no cut, and behaves as a polynomial at infinity (as does uj{Q), i.e. it is a polynomial. 
We thus have the algebraic equation for uj{z): 

S(C) = RiOi e-^ r.,(-^^«i(^) (5,4) 

where T^ is the order r Chebychev polynomial: Tr(cos0) = cosr0. 

Let us now examine the degree of this equation. It is polynomial in u and (, but 
we have to keep in mind that the resolvent of our problem is not uj{(), but W{z) = 
f{z) + gz{z - u), with C = Vfo - f and ^ = 2; - -7j-;^{f + (2 - q)cu). 



c(4-q) ' 



The right-hand side of Eq. (|5.4|) , seen as a polynomial in ( and z, is even and of 
order 2r in (, thus it is also a polynomial in /. 

As for the left-hand side of the equation, S{() verifies S{—() = (—1)''+' S{(), thus, 
when r + / is odd, we have to take the square of Eq. ( |5.4| ) to have our final algebraic 
equation for / which should be of degree 2r. When r + / is even, however, Eq. (|5.4| ) is 
already polynomial in /, and the degree of the equation should be r. 

Moreover, it is easy to check that the order 2r terms on both sides of Eq. ( ^.4|) are 
the same. Finally, the degree d in f{z) of the final equation is: 



d = 2r — 1 if r + I is odd 

d = r — 1 ifr + /is even (5.5) 

while it is of degree d + 1 in z. 

This formula generalizes the presently known results ||^, |^: 

for z/ = 1/3, q = 1, and d = 2 

for u = 1/2, q = 2, and rf = 3 (5.6) 

for z/ = 2/3, g = 3, and d = 5 

In [^, the author investigated the dilute Potts-1,2,3, and 4 cases by the saddle point 
method. These models coincide, when there is no dilution, with the particular case of 
our models with no branching interactions. This gave rise to algebraic equations when 
g 7^ 4 which correspond to our results. However, no general result was given. In M, the 



author used the equations of motion to solve Potts-3 with branching interactions, as 
well as for the resolution of Potts-cxo. But, there again, there was no general expression. 

We are now going to derive the general Potts-g exponent from Eq. (|4.1| ). If, when 
uj{() is singular (i.e. ( is close to one of the bounds of the physical cut of u), uj{—() is 
not, then uj cannot have more complex singularities than half- integer exponents. Thus, 
whereas in the generic case the physical cut [a, b] verifies ab > 0, when the model is 
at the Potts critical point, a (or b) is equal to zero. This means that the bound of 
the unphysical semi-infinite cut we have suppressed by changing variables from / to 
C = V/o — / coincides in that case with the physical cut. 

Let us express: u!{() = C{~()" + regular part Eq. ( [4.1| ) shows immediately that 
the regular part is equal to zero, and that 



^2ina +1 + 2 cos(z/7r) e*™ = i.e. 2-2 cos(z/7r) = g = 2 + 2cos(a7r) 
thus 

, , , , ±v+l+2p 

and 

f r^ {z — const) ii'+i+ap 



As we expect the exponent for / to be greater than one, we have the exponent -^ 
for / and the string exponent, using [|l^]'s formula, is: 

2 (liz^) 



+1 



7. = 1 



±z/ + l (l-±z^) 

6 Conclusion 

In this article, we have obtained general loop equations for the Potts-g model extended 
to branched polymers, which allow us to relate the even and odd parts of the resolvent. 
This relation is then equivalent to an 0(n)-like equation, from which we derived the 
Potts critical exponents and the degree of the algebraic equation which appears in the 
rational case. This last result generalizes the known results for g = 1, 2 and 3 |^, |^ to 
any g = 2 — 2 cos(//r tt) with / < r integers. Moreover, such loop equations as we have 
obtained are necessary when one wishes to apply the renormalization group techniques 
to Potts-g. The study of the renormalization group flows near g = 4 may then provide 
us with useful information about the C = 1 transition. Since g = oo is a real C = oo 
model |T5| and the orders of the phase transitions differ between the fiat and random 



surface Potts models, we may find for large g a set of real C > 1 models. 

Let us also stress that the equations of motion method may be used when A^ is finite 
i.e. on non-planar surfaces, in contrast to the saddle point method. Furthermore, they 
are also less dependent on the analytic structure of the resolvent than the saddle point 
method of p, |^. Finally, as we know how to solve the 0{n) model exactly for general 
values of n, we hope to be soon able to obtain the general expression of the Potts-g 
free energy and operators for any value of g. 
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